Current reversal in collective ratchets induced by lattice instability 
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A collective mechanism for current reversal in superconducting vortex ratchets is proposed. The 
mechanism is based on a two-dimensional instability of the ground state (T = 0) of the system. 
We illustrate our results with numerical simulations and experiments in Nb superconducting films 
fabricated on top of Si substrates with artificially induced asymmetric pinning centers. 
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I. INTRODUCTION 



a substrate with a one-dimensional modulation. 
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Rectification of motion and fluctuations in the 
nanoscale is becoming a major field of research^. Rec- 
tifying mechanisms or ratchets have been used to ex- 
plain how protein motors work^* 3 - and to design new 
separation techniques 4 or synthetical chemical motors 5 -" 6 -. 
Superconductors have become a powerful tool to study 
ratchet mechanisms^. Recently, a superconducting vor- 
tex ratchet device has been reported by Villegas et al£. In 
that experiment, the rocking ratchet mechanism is due to 
a superconducting film patterned with a lattice of asym- 
metric potentials acting as pinning centers. An input ac 
current yields an output dc voltage in the superconduct- 
ing film, revealing a rich phenomenology of single and 
multiple current reversal a 8 ' 9 ' 10 , very sensitive to the un- 
derlying vortex dynamics. The interest of vortex ratch- 
ets is then twofold: it reveals new collective rectification 
mechanisms and sheds light on the physics of vortices in 
superconductors. 

Current reversals in vortex ratchets have been ex- 
plained, in the framework of one-dimensional models, by 
the coexistence of pinned and interstitial vortices mov- 
ing in opposite directions^ or by the interaction between 
vortices within the pinning centers in one-dimensional 
channels 9 . Olson and Reichhardfeii have studied numeri- 
cally a two-dimensional model, obtaining current reversal 
only when interstitial vortices are present in the ground 
state. They provide intuitive explanations of different 
rectification mechanisms based on local interactions be- 
tween vortices and pinning centers. 

In this Letter we experimentally show that current re- 
versal can also occur when interstitial vortices are absent 
in the ground state (T=0). Remarkably, current rever- 
sal disappears increasing either the pinning strength or 
the temperature. We also present numerical simulations 
of vortices as interacting Brownian particles in two di- 
mensions, indicating that this current reversal is due to 
a new collective effect: an instability of the ground state, 
selective to the sign of the applied force. The influence of 
lattice instabilities on rectification has been also recently 
analyzed by Lu et alf 1 ^, for two dimensional vortices in 



II. EXPERIMENTAL RESULTS 

For our experiments, two types of asymmetric pinning 
centers have been fabricated: magnetic (Ni) nanotrian- 
gles and non-magnetic (Cu) nanotriangles. The vortex 
pinning force is enhanced by magnetic centers in com- 
parison with non-magnetic centers 1 ^. The nanotriangles 
were fabricated using e-beam lithography techniques and 
Si (100) wafers as substrates. The Ni or Cu arrays of 
nanotriangles, on top of the substrate, are covered with 
a sputtered Nb thin film. Ni or Cu thickness (triangles 
height) is 40 nm and Nb film is 100 nm thick. Fur- 
ther details on this fabrication technique can be found 
elsewhere 1 ^. For the present work, we have fabricated 
arrays with the same nanotriangle dimensions and array 
periodicity than those in Ref. §• 

We have measured magnetotransport in these films us- 
ing a commercial He cryostat. The variable tempera- 
ture insert allows controlling temperature with stability 
of 1 mK. For these experiments, samples were patterned 
with a cross-shaped measuring bridge^, by using optical 
lithography and ion-etching. This patterned bridge al- 
lows us to control the Lorentz force on the vortices in the 
mixed state: taking into account Ft = J x n(f>Q (with 
(f>o = 2.07 x 10~ 15 Wb and n a unitary vector parallel to 
the applied magnetic field). On the other hand, from the 
expression for the electric field E = B x v, where B is 
the applied magnetic field and v the vortex-lattice veloc- 
ity, we can calculate this velocity v — V/ (dB) from the 
measured voltage drops V (d being the distance between 
contacts). See Ref. 8 for more experimental details. 

The dc magnetoresistance in the mixed state of sam- 
ples with periodic arrays of pinning centers exhibits well- 
known commensurability phenomen a 15 ' 16 , in which min- 
ima develop as a consequence of geometrical matching 
between the vortex-lattice and the underlying periodic 
structure. These minima are equally spaced (two neigh- 
bor minima are always separated by the same magnetic 
field value). For example, in the case of square arrays 
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of nanostructured pinning centers, minima appear at ap- 
plied magnetic fields H m — n(0 o /o 2 )j where a is the lat- 
tice parameter of the square array. Hence, the number of 
vortices n per array unit cell can be known by simple in- 
spection of the dc magnetoresistance R(H) curves. More- 
over, for non-magnetic pinning centers Mkrtchyan and 
Shmidfei£ have given a rough estimation of the maximum 
number of vortices that can be pinned in each center, 
which could confirm the matching field minima values. 
This filling factor can be calculated as the ratio between 
half the dimension of the pinning center (the triangle side 
is around 650 nm)£ and two times the superconduct- 
ing coherence length (around 60 nm for these samples 
and temperatures close to the critical temperature)^. 
In our samples this rough estimation gives us approxi- 
mately three vortices per triangle, in agreement with the 
matching field data (see also |8j) . Therefore, we know, for 
selected values of the applied magnetic field, how many 
vortices there are per unit cell and where they are, this is, 
if they are interstitial vortices or vortices in the pinning 
centers. 



We want to underline that collective behavior of vor- 
tices in films with periodic pinning is crucial to under- 
stand vortex lattice reconfiguration effects 18 or vortex 
channeling^. Also, collective effects in ratchets have 
been shown to yield new interesting phenomena 20 . 



Measurements of the vortex lattice average velocity as 
a function of the applied force are depicted in Fig. QJa) 
for sample (A), Ni triangles, and three vortex per triangle 
(n = 3). There is no current reversal in the rectified 
signal in agreement with the explanation given in Ref. 8 
(see below). Experiments for sample (B), Cu triangles, 
also for n = 3 (Fig. HJb)) show a similar behavior for 
T = 8.165K ( T = 0.99 T c ), but a current reversal appears 
decreasing the temperature, despite the fact that there 
are still three vortices per triangle. 



In the one-dimensional approach reported in Ref. [H, 
the vortex lattice does not play any role, and the current 
reversal was explained assuming that pinned and inter- 
stitial vortices are rectified in opposite directions. How- 
ever, this assumption fails to explain current reversal for 
the matching field n — 3, since in this case there are 
no interstitial vortices in the ground state of the system. 
The simulations performed by Olson and Reichhardt 11 
do not exhibit any current reversal for n = 3 and the 
rectification mechanisms that they propose do not apply 
for this case either. Therefore, a new explanation of cur- 
rent reversal is required. We have found by numerical 
simulations such an explanation based on the interplay 
between the vortex lattice and the geometry of the tri- 
angular defects. 
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FIG. 1: Net velocity of vortices versus ac Lorentz force am- 
plitude (ui = 10kHz). Array periodicity 770 nm, triangle 
base 620 nm. (a) Sample (A): Nb film with Ni triangles 
T c = 8.35iC~, T/T c = 0.99 (T = 8.265K) (□), T/T c = 0.98 
(T = 8.174K) (o). (b) Sample (B): Nb film with Cu triangles 
T c = 8.24K. T = 8.165K (□), T = 8.102K (x), T = 8.08K 
(o), T = 8.065K (0). The negative velocity part of each curve 
is shown in the inset. 



III. THEORETICAL MODEL 

A. Simulations 

The simulations have been performed by numerically 
solving Langevin equations for the movement of the vor- 
tices 

V ±= ~d x U vv + VV p + F cxt + £{t). (1) 

F ext is the Lorentz force resulting from the applied 
current, U vv the usual vortex- vortex interaction (see 
Refill), V p the pinning potential, r\ the friction coeffi- 
cient, and £(t) a white gaussian thermal noise. The pin- 
ning force and vortex- vortex interaction must be in agree- 
ment with the experimental situation of three vortices per 
pinning site (triangle). The interaction of the vortices 
with the pinning defects is modelled by a potential V p 
in the shape of a triangle and with a hyperbolic tangent 
profile. The depth of the potential is V p q = 0.002 pN/im, 
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FIG. 2: Simulation results. Three vortices per triangle, pin- 
ning 0.002 pN /jm. Inset: pinning 0.0036 pN/xm. 



so that depinning forces have typical values of the order 
of 10~ 6 N/m. Experiment s 2 ^ 23 show that the system is 
adiabatic in the region of frequencies used. This allows 
us to obtain the expected ac signal from the velocity- 
force response curves obtained in simulations, both for 
constant positive and negative applied force. 

Results are shown in Fig. [5] where a window of down- 
ward rectification can be observed. Let us stress at 
this point that the "natural" direction of rectification 
for pinned vortices is upwards. The reason is that the 
vortex feels a smaller force at the tip of the triangle than 
at the triangle base. The force at the tip can be, for an 
equilateral triangle, as lower as half the force at the base 
of the triangle. 



B. Ground state instability 

What is then the origin of the downward rectification 
for small forces? Our simulations indicate that the rec- 
tification is due to a two-dimensional instability of the 
ground state under small downward forces. The ground 
state of the system consists of vortices located in the cor- 
ners of each triangle, without interstitial vortices. How- 
ever, for finite temperature, there are some "defective tri- 
angles" in the actual configuration of vortices, i.e., some 
triangles which have only two vortices inside. Accord- 
ingly, some interstitial vortices will be randomly spread 
along the sample. Fig. [3] shows one of these configura- 
tion with only one defective triangle out of 6 x 6. For low 
forces, motion is induced by the interstitial vortices both 
in the downward and upward direction. However, con- 
trary to the picture presented in previous works^, there is 
not a continuous motion of interstitial vortices along the 
space between triangles: the interstitial vortex enters the 
nearest triangle expelling one of the three vortices inside. 

The aforementioned instability, selective to the sign of 
the external force, is shown in Fig. [3] We have chosen an 
initial condition with only one interstitial vortex and one 
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FIG. 3: Downward rectification mechanism at T = K: snap- 
shots from simulations. Initial condition (left) and configu- 
rations after evolution (r = 6.5 x 10~ s) with positive force 
F — 0.032pN (right, up) and negative force — F (right, down). 
Dashed circles show the defective triangle and the only inter- 
stitial present in the initial condition. Thick arrows point out 
the only two columns presenting any motion when positive 
force is applied. Comparison of the right panels shows that 
the number of interstitial vortices ( x ) dramatically increases 
when the external force points downwards. 



defective triangle, located at some distance. In the right 
figures we have plotted the configurations of vortices af- 
ter an upward and downward force has been applied for 
t = 6.5 x 10~ 9 seconds (long enough for a depinned vor- 
tex to cover the whole sample several times), respectively. 
This simulations evolved at zero temperature to show the 
mechanism in a clearer manner. We see that, in the case 
of the upward force, the interstitial vortex remains in one 
single column. The column with a defective triangle also 
moves, but again vortices remain in that column. As a 
consequence, there is a positive current of vortices but 
the motion is constrained to two columns. When the in- 
terstitial vortex enters a triangle, the top vortex in the 
triangle moves upwards, out of it and into the following 
one. In addition to that, in the defective triangle, one 
of the two vortices always escapes through the tip, be- 
coming and interstitial vortex and triggering a process 
similar to that in the column with an extra vortex. Both 
motions propagate along the column without disturbing 
the neighboring columns. 

On the other hand, under a downward force, the ini- 
tial defect in the ground state propagates along the whole 
sample, as it is clearly shown in Fig. [3] A more detailed 
analysis of the simulations show that, when an interstitial 
vortex enters a triangle, one of the two bottom vortices 
is expelled but now can move to a triangle in one of the 
nearest columns. It even happens frequently that the two 
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FIG. 4: Vortex-vortex interaction versus vortex-vortex dis- 
tance. Solid line: interaction used in simulations corre- 
sponding to current reversal (see Fig. [2] pinning V p o = 
0.002pN/im.). Dashed line: interaction used in simulations 
showing no current reversal, which are depicted in the inset 
(pinning V p0 = 0.002pN /Ltm). 

vortices in the base of the triangle are depinned. Con- 
sequently, the initial defect is then spread out along the 
horizontal direction, yielding a considerable large fraction 
of depinned vortices which increase the overall motion in 
the system, yielding a net downward rectification and the 
corresponding current reversal. 

As was noted before, current inversion may disappear 
when the pinning potential is increased. This effect was 
observed in experiments where the Cu triangles were re- 
placed by Ni ones which have a higher pinning strength 
(see Fig. QJa)), arL d it is also reproduced in our simula- 
tions. Results for pinning potential V p o = 0.0036pN^m 
do not show current reversal, as depicted in the inset 
in Fig. O The increased pinning strength implies that 
a larger force is needed to depin vortices when pushing 
downwards. According to simulations, when the same 



force is applied upwards, the aforementioned columnar 
movement starts, but also neighboring columns without 
defective triangles or extra vortices depin, and eventually 
motion spreads all over the sample. In other words, for 
these moderate values of the force, the instability appears 
in both directions. 

Finally, current reversal in sample (B), Cu triangles, 
also vanishes when temperature is raised (see Fig.QJb)). 
We can explain this behavior taking into account that the 
penetration depth of the superconductor increases with 
temperature. As a result, the vortex-vortex interaction 
strength, U vv , decreases at short distance but its range 
becomes longer—, as depicted in Fig. [4] In this case we 
observe that the long range vortex lattice order precludes 
the instability responsible for current reversal (see the 
inset in Fig. [4|). 

IV. CONCLUSIONS 

In summary, our simulations show that the asym- 
metric substrate induces an instability sensitive to the 
direction of the external force, affecting rectification. 
This effect can explain the current reversal observed in 
our experiments for n — 3. Moreover, our work indi- 
cates that the lattice configuration, and consequently 
its dynamical properties, can be controlled by exter- 
nal forces. This interplay between rectification, driving 
forces and lattice configuration can induce other inter- 
esting phenomena such as transitions between different 
lattice configurations 1 -, and could help to design new 
rectifying devices, not only in superconducting films, but 
also in other two dimensional collective systems, such as 
Josephson arrays or colloidal suspensions. 
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